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Abstract. The application of the virial theorem provides a tool to estimate supermassive 
black hole (BH) masses in large samples of active galactic nuclei (AGN) with broad emis- 
sion lines at all redshifts and luminosities, if the broad line region (BLR) is gravitationally 
bound. In this paper we discuss the importance of radiation forces on BLR clouds arising 
from the deposition of momentum by ionizing photons. Such radiation forces counteract 
gravitational ones and, if not taken into account, BH masses can be severely underestimated. 
We provide virial relations corrected for the effect of radiation pressure and we discuss their 
physical meaning and application. If these corrections to virial masses, calibrated with low 
luminosity objects, are extrapolated to high luminosities then the BLRs of most quasars 
might be gravitationally unbound. The importance of radiation forces in high luminosity 
objects must be thoroughly investigated to assess the reliability of quasar BH masses. 
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1. The BH mass ladder 

The discoveries of the tight correlations be- 
tween the mass of BHs and luminosity, mass, 
velocity dispersion and surface brightness pro- 
file of their host spheroids have produced an 
enormous impact. MeH-host galaxy relations 
allow a "census" of BHs in the local universe 
showing that local BHs were grown during 
AGN activity. Most of all, M B H-host galaxy 
relations indicate that BHs are an essential 
element in the evolution of galaxies and the 
link between BH and host galaxy is probably 
caused by AGN feedback. It is fundamental to 
measure BH masses in large samples of ob- 



jects from zero to high redshifts to fully un- 
derstand the effects of BH growth on galaxy 
evolution. Direct BH mass estimates based on 
spatially resolved stellar and gas kinematics 
are possible only in the local universe but 
can be the starting point to calibrate less di- 
rect techn i ques following a "BH mass ladder" 
dPetersonl |2004|) . Rung one of the BH mass 
ladder is provided by gas and stellar kine- 
matical measurements resultin g in the Mrh-cTr 
and M fiH-^sph relations (e.g. iFerrarese & Fordl 
2005). Rung two is provided by reverberation 
mapping (RM) measurements; BH masses are 
given by M BH = /AV 2 /? B lr/G, where R BLR 
is the BLR average distance from the BH es- 
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timated with reverberation mapping, AV is the 
width of the broad emission line and / is a scal- 
ing factor which depends on (unknown) BLR 
properties. Although this technique is poten- 
tially plagued by many unknown systematic er- 
rors, BH masses from RM are in agreement 
with the MBH-c e relation of normal galax- 
ies. This second ru ng has been calibrated by 
lOnken et al.l (120041) who estimated / assum- 
ing that the AGN with RM observations fol- 
low the Meti-Ce relation. In principle this tech- 
nique is not limited by distance but, in prac- 
tice, RM observations of high-z, high-L AGNs 
are exceedingly demanding in ter m of observ- 
ing ti me. The ^blr - L relation (Kaspi et al., 
2000) shows that continuum luminosity acts as 
a proxy for /?blr leading to the so-called single 
epoch (SE) virial mass estimates. Rung three of 
the BH m ass ladder has been calibrated, a mong 
others, by Vestergaar d & Peterson! f2006) start- 
ing from RM measurements. Since a single 
measurement can be wrong by up to a factor 
of ~ 10, SE Mbh are considered accurate only 
in a statistical sense and are routinely used to 
estimate Mbh in large sample of AGNs from 
zero to high redshifts. 



2. The effect of radiation pressure 

One of the basic assumptions of reverbera- 
tion mapping is that the B LR is photoion- 
ized dBlandford & McKeel Il982l) implying 
that BLR clouds are subject to radiation forces 
arising from ionizing photon momentum depo- 
sition. Therefore BLR clouds will effectively 
"see" a smaller BH mass. An order of magni- 
tude estimate of this radiation force can be ob- 
tained with a simple heuristic model in which 
BLR clouds are test particles, optically thick 
to ionizing photons. The corrected "virial rela- 
tion" is then 



fV z = - M BH -— — 



(1) 



where r is the cloud distance from the 
AGN/BH, TO p is the proton mass, a = Lj on /L 
is the ionizing to total luminosity ratio and Nu 
is the total cloud column density in the direc- 
tion of the AGN/BH (see Mar coni et aDl2008l 



hereafter M08, for more details). The corrected 
mass estimator is 



V 2 r 
Mm = f^r- + 



G 4jtGc m p Nh 



(2) 



the first term is the "canonical" virial relation 
while the second one is the radiation force 
correction, which can easily become relevant. 
For example, in AGNs with moderate column 
densities (Nu — 10 23 ) and luminosities (L 
10 11 L Q ), the correction becomes of the order 
of 2.7 x 10 7 M , a typical value of the BH mass 
in a Seyfert galaxy (we have assumed a = 0.6). 
The above correction to Mbh is model depen- 
dent and Nu is unknown but we can consider a 
virial estimator of the form: 



Mi 



BH 



V 2 r 
f— + 8 



AL A ( 5 100 A) 



10 44 ergs- ! 



(3) 



One can then repeat 
rungs 2 and 3 as in 



the calibrati o ns for 
Onken et all (2004) 



and lVestereaard & Peterson! d200 6). The whole 
processes is presented and discussed in M08 
and here we only outline the main results. 
Calibration of Rung 2 is performed by impos- 
ing that the galaxies with RM data fall on the 
MBH-0" e relation and the virial estimator for re- 
verberation mapping observations (RM) is: 



Mr 



= (3.1 ± 1.4)^ + (l0 7 - 6±0 - 3 L 510 o)M £ 



where V 2 R/G is the virial product (in units 
of M»: iPeterson et al.ll2004h and L 5 i o is AL A 
at 5100 A (in units of 10 44 ergs). The value 
/ = 3.1 here is alm ost a factor 2 smal ler than 
/ = 5.5 found by lOnken et all (120041) due to 
the radiation pressure correction. The ~ 0.5 dex 
r.m.s. scatter of BH masses around MBH-0~ e 
provides an indication of systematic errors af- 
fecting Mbh(RM) estimates compared to gas 
or stellar kinematical measurements. The in- 
troduction of the g parameter is not made to 
improve the fit (the scatter is not significantly 
reduced) but is required by the physics of BLR 
clouds. Our ability to determine an accurate 
empirical value of g is limited, as were previ- 
ous efforts to determine /, by size, composi- 
tion and accuracy of the existing RM database. 
In particular it contains few sources with high 
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Eddington ratios, which provide the tightest 
constraints on g. With this caveat in mind, 
however, we can jump onto the next rung of 
the BH mass ladder. Calibration of Rung 3 is 
performed by imposing that SE Mbh estimates 
in galaxies with RM data are in agreement with 
the corresponding Mbh(RM) and the virial es- 
timator from single epoch (SE) observations is: 

M BH /M = 10" 3 S^ L o.5 o + 10 "2!™ L5100 

where V H/j is the FWHM of broad Hft (in 
units of 1000 km s). Without the radiation pres- 
sure correction, M B h(SE)/M B h(RM), has a 
r.m.s. scatter of ~ 0.4 dex. Significantly, the 
scatter drops to ~ 0.2 dex when radiation pres- 
sure is taken into account indicating that SE 
masses are more accurate M B h estimators than 
previously thought (see Fig. 2 of M08). This 
0.2 dex r.m.s. scatter indicates only system- 
atic errors made when using SE instead of RM 
masses. The physical meaning of this empiri- 
cal g value (g ^ 10 77 M ) can be assessed with 
our heuristic model. From Eq. [2] the average 
Nn needed to obtain the empirical g value is 
TYh - 1.1 X 10 23 cirT 2 and is remarkably similar 
with the results from photoionization model- 
ing studies of the BLR. Our empirical determi- 
nation has indeed a simple physical meaning 
consistent with our current knowledge of the 
BLR. A further indication of the validity of our 
corrected virial relation is provided by Narrow 
Line Seyfert 1 galaxies. These galaxies are be- 
lieved to have small BHs compared to other 
AGNs an d what is expected from the Mnu-<r„ 
relation dGrupe & Mathurl l2004t IZhou et all 
2006). Their small BH masses and high L/L E dd 
ratios suggest that these galaxies are rapidly 
growing their BHs. With the radiation pressure 
correction, NLSls move on to the MeH-fe re- 
lation, indicating that their BH masses were 
likely underestimated (see M08 for more de- 
tails). 

3. High luminosity objects 

Virial relations involving UV lines are im- 
portant to study sources at high redshift and 
it is possible to obtain virial relations for 
Civ (A 1549 A) and Mgn (,12798 A) which 
take into account radiation pressure like for 



Hy6 (Marconi et al. 2008, in preparation). 
Preliminary relations are: 

M BH /M - 10 6 - 5 V 2 IV L? 3 5 50 + 10 7 - 2 L 1350 

M BH /M - 10 6I y^ gU L° 5 00 + 10 7 - 5 L 30 oo 

where V C iv, Vivigii are the FWHM of Civ 
and Mg n broad lines respectively (in units of 
1000 kms~') and Z4350, £3000 are the continua 
at 1350 and 3000 A respectively (same no- 
tation as £5100)- The differences in g factors 
w.r.t. Hy6 depends on the different bolomet- 
ric corrections. To study the effect of radia- 
tion pressure corrected virial relations on high 
lu minosity ob j ects w e consider the catalogue 
bv lShen et al.l (f2008) which contains line mea- 
surements for ~ 60000 SDSS quasars. In Fig.Q] 
we plot the ratio of new-to-old BH masses as 
a function with redshift (left) and the distribu- 
tion of L/Lem ratios as a function of luminos- 
ity (right). These two plots outline a few im- 
portant points. BH masses corrected for radi- 
ation pressure are not necessarily larger than 
older ones; the average ratio ranges from ~ 1 
to ~ 2 at high z with a broad distribution of 
values. As shown above, the / factor computed 
considering the radiation pressure correction is 
smaller by almost a factor 2; objects with low 
luminosities compared to their virial products 
will have a negligible radiation pressure cor- 
rection and a BH mass which can be a factor 
~ 2 smaller. At high z, sources in the catalogue 
have larger average luminosities while the dis- 
tribution of li ne widths does not change appre- 
ciably with L dShen et al.Ll2008l) : hence the ra- 
diation pressure correction is more important 
and BH masses are, one average, larger. The 
distribution of L/Z-Edd with luminosity is pecu- 
liar, reaching a well defined limiting value at 
about L/^Edd - 0.15. Indeed when the lumi- 
nosity is large the virial product becomes neg- 
ligible with respect to the radiation correction. 
Hence Mbh k L and the L/^Edd value saturates 
to a constant L*/LEdd value, where L* is the 
critical luminosity which unbinds BLR clouds. 
This behaviour is a clear consequence of the 
adopted correction to the virial relation (Eq.[TJi. 
Moreover, such a relation is correct only if the 
right side is positive, i.e. if the luminosity is 
smaller than and the BLR is gravitationally 
bound, the fundamental underlying assumption 
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Fig. 1. Left: ratio between virial BH masses corrected (NEW) and non-corrected (OLD) for radiation pres- 
sure. Red filled squares denote average values in given redshift bins. Center: luminosity distribution of 
L/L Edd from NEW virial masses. Right: Contours of constant line FWHM (solid black lines) in the M BH - L 
plane. Numbers denote FWHM values in km s . The red and green lines denote the loci for L = L^,n and 
L = L t , where L t is the critical luminosity at which the BLR is gravitationally unbound. 



which, of course, cannot be tested a-posteriori 
using the derived Mbh value. Fig. QJ) clearly 
indicates that many high L sources are close to 
the critical luminosity and this raises a poten- 
tial problem. If at quasar luminosities the BLR 
is experiencing only a small fraction, e.g. < 
10%, of the real gravitational potential then 
one should wonder whether the assumption of 
a gravitationally bound BLR is still acceptable. 
In Fig. [TJ; we consider the Mbh - L plane. For 
each point in this plane one can estimate the ex- 
pected line width using Eq.[T]and the Rblr - L 
relation. Black lines denote constant FWHM 
values while the red and green lines represent 
the L = Lndd (solid) and L — L* loci. The area 
between the red and green lines is populated by 
sources which are below Eddington but have 
their BLRs unbound. Since most observed 
quasar spectra hav e FWHM < 10000 km/ s 
dShen et all 120081) . a quasar with L > 10 12 L o 
will always have BLR clouds close to be- 
ing unbound. A similar situation characterizes 
NLS 1 at lower luminosities ( 1 000 < FWHM < 
2000 km/ s) which, in this picture, are low lu- 
minosity analogues of high luminosity quasars 
in terms of their BLR. The situation is differ- 
ent for Seyfert galaxies which populate a re- 
gion where BLR clouds are well bound grav- 
itationally. This simple graph outlines a po- 
tential problem: the radiation force on BLR 
clouds in high L quasars is extremely impor- 
tant compared to the gravitational attraction 



posing a potential problem on the validity of 
virial relations. Of course, BLR clouds in high 
L quasars could be dominated by gravity if for 
instance the continuum is not isotropical and 
the BLR is i lluminated by a m uch fainter con- 
tinuum (e.g. Proga et al. 2008), or BLR clouds 
in quasars might have much larger column den- 
sities (e.g. Nh > 10 25 cm~ 2 ) than in lower lumi- 
nosity objects. Overall, it is clear that the radi- 
ation force is an important effect which should 
be taken into account in virial masses and that 
it is mandatory to assess whether the virial as- 
sumption of a gravitationally bound BLR is 
still valid in high luminosity quasars. 
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